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1 Introduction 

In paper [1] E.J. Parkes presented a categorization of solutions of the equation dubbed the 
extended reduced Ostrovsky equation (exROE). The equation studied has the form 



0, 



where T> 



d_ d_ 
dt dx 



(1) 



with and (3 being constant coefficients. This equation was derived from the Hirota- 

Satsuma-type shallow water wave equation considered in [2] (for details see [1]). 

For stationary solutions, i.e. solutions in the form of travelling waves depending only on one 
variable x = ^ — — ^o, this equation reduces to the simple third-order ODE: 
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(2) 



where V stands for the wave speed and w = u — V. (Note, that in many contemporary papers 
including [1] authors call such solutions simply "travelling- wave solutions". Such terminology 
seems not good as nonstationary propagating waves also are travelling waves. The term "sta- 
tionary waves" widely used earlier seems more adequate for the waves considered here.) In 
paper [1], equation (2) was reduced by means of a series of transformations of dependent and 
independent variables to an auxiliary equation whose solutions were actually categorized subject 
to some restrictions on the equation coefficients, viz.: 



p + q^O, qV-/3^0 



(one more restriction on the constant of integration for that auxiliary equation, S = 0, was 
used in [1]). Under these restrictions, solutions to equation (2) were found in analytical form 
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and corresponding wave profiles were illustrated graphically. Among solutions obtained there 
are both periodic and solitary type solutions including multivalued loop periodic waves and 
loop-solitons. 

Similar loop solutions to exROE and some other equations were earlier obtained by Ji-Bin 
Li [3] who came to the conclusion that loop solutions actually are compound solutions which 
consist of three different independent branches. These branches may be used in various combi- 
nations representing several types of stationary propagating singular waves (waves with infinite 
gradients). This conclusion completely coincides with the conclusion of paper [4] where a com- 
plete classification of stationary solutions of ROE was presented. ROE derived by L.A. Ostrov- 
sky [5] in 1978 as a model for the description of long waves in a rotating ocean (see [4] and 
references therein) can be treated as a particular case of exROE with p = q and = (see [1]). 

Below an analysis of stationary solutions to equation (2) is presented by the direct method 
avoiding any redundant transformations of variables. The method used is based on the phase 
plane concept and analogy of the equation studied with the Newtonian equation for the point 
particle in a potential field. Such approach seems more vivid and free of aforementioned re- 
strictions. This work can be considered also as complementary to paper [1] as the analysis 
presented may be helpful in the understanding of basic properties of stationary solutions of 
equation (2). 



2 Mechanical analogy, potential function 
and phase-plane method 

Equation (2) can be integrated once resulting in 

+ ^{p + q)w^ + {pV + f3)w = Ci, (3) 

where Ci is a constant of integration. By multiplying this equation by dw/dx and integrating 
once again, the equation can be reduced to the form of energy conservation for a point particle 
of unit mass moving in the potential field P{w)\ 

1 / dw \ ^ 



where the effective "potential energy" as a function of "displacement" w is 



6 w w"^ 



and C2 is another constant of integration. The constant E = —{pV + f3)/2 plays a role of 
the total energy of the particle, i.e. the sum of the "kinetic energy", K = {l/2){dw / dxY ^ and 
the "potential energy", P{w). As follows from equation (4), real solutions can exist only for 
E > P(w). Various cases of the potential function (5) are considered below and corresponding 
bounded solutions are constructed. Unbounded solutions are not considered in this paper as 
they are less interesting from the physical point of view; nevertheless, their qualitative behavior 
becomes clear from the general view of corresponding phase portraits. 



3 Particular case: p -\- q = 

Consider first a particular case when the coefficients in equation (1) are such that p + g = 0. 
Note, this is one of the cases which were omitted from the consideration in paper [1]. The 
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Figure 1. a) Potential function for the case p + g' = 0, C2 = and two values of Ci: Ci = 1 (solid lines), 
and Ci = — 1 (dashed lines). Two dots illustrate possible motions of point particles in the potential field, 
b) Phase plane corresponding to the potential function with Ci = 1 and different values of E. Line 1: 
E = -1; line 2: E = -0.5; lines 3: E = 0.5; lines 4: E = 1. 



potential function (5) simplifies in this case. However, a variety of subcases can be distinguished 
nevertheless even in this case depending on the coefficients Ci and €2- All these subcases are 
studied in detail below. 

3a. If C2 = 0, the potential function represents a set of antisymmetric hyperbolas located 
either in the first and third quadrants or in the second and fourth quadrants as shown in Fig. la. 
The corresponding phase plane {w^w')^ where w' — dw/dx, is shown in Fig. lb for Ci — 1 (for 
Ci — —1 the phase plane is mirror symmetrical with respect to the vertical axis). For other 
values of Ci phase portraits are qualitatively similar to that shown in Fig. 1 for Ci = 1. 

Analysis of the phase portrait shows that there are no bounded solutions for any positive £"; 
corresponding trajectories both in the left half and right half of the phase plane go to infinity 
on w (see, e.g., lines 3 and 4 in Fig. lb). Meanwhile, solutions bounded on w do exist for negative 
values of E (i.e. for V > —f^/p), but they possess infinite derivatives when w = 0. Consider, 
for instance, motion of an affix along the line 2 in Fig. lb (Ci = 1) from = 00 towards the 
axis w where = 0. The qualitative character of the motion becomes clear if we interpret it in 
terms of "particle coordinate" w and "particle velocity" treating ^ as the time. The motion 
originates at some "time" ^0 with infinite derivative and zero "particle coordinate" w = 0. Then, 
the "particle coordinate" w increases to some maximum value u^max — 

-Ci/E {E < 0) as the 

"particle velocity" is positive. Eventually it comes to the rest having zero derivative = 
and w = tt;max- Another independent branch of solution for the same value of E corresponds 
to the affix motion along the line 1 from the previously described rest point at axis w towards 
= — oc and w = 0. 

All bounded analytical solutions for this case can be presented in the universal implicit form: 



arctan (^^^3^) ~ VvC^-v) 



(6) 



where y = -Ew/Ci, ^ = -V2{-Ey/^x/Ci and is an arbitrary constant of integration. This 
solution consists of two independent branches which correspond to signs plus or minus in front 
of the square brackets in equation (6). Each branch is defined only on a compact support of 
axis ^: either on —7t/2 < ^ — < or on < ^ — ^ (see lines 1 and V in Fig. 2). 
With the appropriate choice of constants one can create a variety of different solutions, e.g., 
the y-shape wave (see lines 2 and 2'), or a smooth-crest compacton, i.e. a compound solitary 
wave defined only for |C — Co| ^ (see lines 3 and 3^). Using a translational invariance of 
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Figure 2. Various particular solutions described by equation (6). 




I / 

■■ -5 ^ 

Figure 3. Maximum of the compacton solution (6) against speed in the original variables, equation (7). 
Dashed vertical line corresponds to the limiting value of V = —/3/p. The plot is generated for Ci = p = 
(3 = 1. 



solutions and their independency of each other, one can create periodic or even chaotic sequences 
of compactons randomly located on axis ^. 

The maximum of the function y(^), I/max = 1, corresponds in terms of w to Wma^ = —Ci/E. 
Using the relationship between w and the original variable u (see above), as well as the definition 
of the constant E, one can deduce the relationship between the wave extreme value (wave 
maximum) and its speed: 

Taking into account that we consider the case of Ci = 1, and negative values of E are possible 
only when V > —jS/p^ the plot of i/max(V^) is such as presented in Fig. 3. 
As follows from equation (7), a wave is entirely negative (i/max < 0), when 



v<l-UW^Spc[-^) 



2p 

provided that p < (8(7i). At greater values of the wave profile contains both positive and 
negative pieces, and for certain value of V the total wave "mass" I — j u(x)dx vanishes (the 
integral here is taken over the entire domain where function u{x) is defined). 
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a) b) 

Figure 4. a) Potential function for the case p + g = 0, Ci = and two values of C2: C2 = 1 (solid lines), 
and C2 = —I (dashed lines). Two dots illustrate possible motions of point particles in the potential field, 
b) Phase plane corresponding to the potential function with C2 = 1 and different values of E. Line 1: 
E = —1; line 2: E = —0.5; lines 3: ^ = 0.5; lines 4: E = 1. All lines are symmetrical with respect to 
axis and are labelled only in the left half of the phase plane. 



3b. A similar analysis can be carried out for the case when Ci = 0, C2 7^ 0. The potential 
function in this case represents a set of symmetric quadratic hyperbolas located either in the first 
and second quadrants or in the third and fourth quadrants as shown in Fig. 4a for C2 = ±1. The 
corresponding phase plane is shown in Fig. 4b for C2 = 1 only (there are no bounded solutions 
for C2 = — 1, therefore this case is not considered here). For other positive values of C2 phase 
portraits are qualitatively similar to that shown in Fig. 4b. 

Analysis of the phase portrait shows that there are no bounded solutions for C2 = — 1, as 
well as for C2 = 1 and any positive E (see, e.g., lines 3 and 4 in Fig. 4b); they exist however for 
C2 = 1 and negative values of E, but possess infinite derivatives at some values of x- iioi"- 
malized variables y = {—E/C2)^^'^w, ^ = —E{2/C2)^^'^X possible solutions can be presented 
in terms of independently chosen function branches describing a unit circle in one of the four 
quadrants, i.e. 

i^-^of + y^-i, (8) 

where ^0 is an arbitrary constant of integration. 

Playing with the constant ^0 one can create again a variety of compacton-type solutions 
including multi- valued solutions. Some examples of solitary compacton solutions are shown in 
Fig. 5a; they include A^-shaped waves, multi-valued circle-shaped waves and semicircle positive- 
polarity pulses (due to symmetry, the polarity of the first and last waves can be inverted). In 
addition to those, various periodic and even chaotic compound waves can be easily constructed; 
one of the possible examples of a periodic solution is shown in Fig. 5b. Each positive or negative 
half-period of any wave consists of two independent branches originating at = and ending 
at 7/ = ±1. The same is true for the pulse-type solutions shown in Fig. 5a; they consist of 
independent symmetrical branches as shown, for example, for the semicircle pulse in Fig. 5a 
where they are labelled by symbols 1 and 2. 

The maximum of the function y{^), ^max = 1, corresponds in terms of w to the wave maxi- 
mum, it;max = {—C2/E)^/'^. Using a relationship between w and the original variable u (see 
above), as well as definition of the constant E^ one can deduce the relationship between the 
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Figure 5. a) Some examples of pulse-type waves described by equation (8): A'^-shaped wave; circle wave 
and semicircle compacton. b) One of the examples of a periodic wave with infinite derivatives at y = 0, 
^ = 2n + 1, where n is an entire number. 



wave maximum and its speed: 



v + 



-C2 



v + 



I 2C2 
pV + 13' 



(9) 



The plot of itmax(^) is presented in Fig. 6 for y > —fi/p in accordance with the chosen constant 
C2 = 1 and £' < 0. 

As follows from equation (9), wave maximum (minimum) cannot be less than the certain 
value, ?7max (— ^^min), which occurs at some speed Vi, where 



P 



,2\ 1/3 



/3 



p 



C2P' 



2\ 1/3 



For all possible values of wave maximum i/max > ^^max, two values of wave speed are pos- 
sible, i.e. two waves of the very same "amplitude" can propagate with different speeds. This is 
illustrated by horizontal dashed line in Fig. 6 drawn for i^max = 2.5. The same is true for waves 
of negative polarity. 

3c. Consider now the case when both Ci and C2 are nonzero but p + g is still zero. There 
are in general four possible combinations of signs of the parameters Ci and C2'. 

i) Ci > 0, C2 > 0; ii) Ci < 0, C2 > 0; iii) Ci > 0, C2 < 0; iv) Ci < 0, C2 < 0. 

The shape of the potential function P{w) and corresponding solutions are different for all 
these cases. However, among them there are only two qualitatively different and independent 
cases, whereas the two others can be obtained from those two cases using simple symmetry 
reasons. This statement is illustrated by Fig. 7, where the potential relief is shown for all four 
aforementioned cases i)-iv). 

As one can see from Fig. 7, cases i) and ii), as well as iii) and iv), are mirror symmetrical 
counterparts of each other with respect to the vertical axis. This implies that solutions for the 
cases i) and ii), and correspondingly, iii) and iv), are related by the simple sign interchange 
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Figure 7. Potential relief for the four different cases, i)-iv), of various signs of constants Ci and C2. 
The plot was generated for Ci = ±1, C2 = ±0.1. 



operation, i.e. W]^^ — — w-^y^^^ — —w-^^y Therefore, below only two qualitatively different cases 
are considered in detail, namely the cases i) and iii). 

Case i) is characterized by an infinite potential well at the origin, w — {). This singularity 
in the potential function corresponds to the existence of a singular straight line it; = on the 
phase plane (see Fig. 8). On both sides from this singular line there are qualitatively similar 
trajectories which correspond to bounded solutions having infinite derivatives at the edges. 
Quantitative difference between the "left-hand side solutions" and "right-hand side solutions", 
apart of their different polarity, is the former solutions (of negative polarity, w < {)) exist for 
E < Pmax, whereas the latter ones (of positive polarity, w >{)) exist for £" < 0. The potential 
function has a maximum 

-Pmax = Cl/{AC2) &tw = -2C2/C1. There are no bounded solutions 

for E > Pmax. 

Consider first bounded solutions which correspond to trajectories shown in the left half- 
plane, It; < 0, in Fig. 8. For a positive value of the parameter E in the range < < Pmax, the 
analytical solution can be presented in the form 



ay) = ±2yg 



i/ iy + 2Qy - 4g(g - 1) - 2g m ^ + + S^T ^ 



(10) 



where ^ = xv^lCi/Ca)^, y = w{Ci/C2), Q = Cf/{AC2E). 
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Figure 8. Phase portrait of equations (3), (4) for the case i) (only those trajectories are shown which 
correspond to particle motion within the potential well in Fig. 7a). Line 1: E = 2.5; line 2: E = 1; 
lines 3: E = -1: line 4: E = -2: line 4: E = -5. 



The range of variability on ^ is: 

lel < 4Q {l - VQln [(V^ + - l] } , 

whereas y varies in the range 

-2 [q - VQ{Q - 1)] <y<0- 

The relationship between the wave minimum and its speed is: 
Ci 



Ur, 



= v + 



+ 



2C2 



pV + /3 V pV + (3 [2C2{pV + (3) 



CI 



+ 1 



(11) 



where pV^ + /? < as > 0. 
If £^ < 0, then the solution is 



\/4Q(Q 



+ 2Q arctan 



(y + 2Q)2 



+ 7rQ 



(12) 



y4Q(Q-l)-(y + 2Q)2^ 

The range of variability on ^ is: |^| < —4(5 [l + \f—Q (arctan \/—Q — 7r/2)] , whereas y varies 
in the range —2 Q + \/ Q{Q — 1) < y ^ 0. The relationship between the wave minimum and 
its speed is also given by equation (11), but with pV + /? > 0. 

Two special cases of solution (10) can be mentioned. When Q = 1 (£" = ^max), solution (10) 
with the appropriate choice of the integration constant reduces to 



m = ±4 



L2 V 2 



(13) 
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This solution is unbounded on ^, i.e. it is defined in the range: |^| < oc. However, the 
solution is bounded on y: —2<y<0. The relationship between the wave minimum and its 
speed is simple as both of them are constant values in this special case: 

Another special case corresponds to Q = oo = 0); in this case equation (10) after appro- 
priate choice of integration constant reduces to: 

2 



e(y) = ±-y^(y-2). (15) 

The range of variability on ^ is: |^| < 4/3, whereas y varies in the range: — 1 < y < 0. The 
relationship between the wave minimum and its speed is also very simple as both of them are 
again constants but different from those given by equation (14); in this case they are: 

p Ci \p Ci 

Bounded solutions corresponding to the trajectories shown in the right half-plane in Fig. 8 
with w > 0, exist only for negative E; they are given by equation (12), but with the different 

range of variability of y: < y < 2 ^—Q + \/Q{Q — 1) . The relationship between the wave 
maximum and its speed is given again by equation (11) where i^max should be substituted instead 
of iXmin and pV + > as £^ < for these solutions. 

Solutions (10), (12), (13) and (15) are shown in Fig. 9. All these solutions are of the com- 
pacton type; they consist of two independent branches which can be matched differently or 
unmatched at all. Lines 2 and 2^ represent an example when two branches are matched so that 
they form a semi-oval; lines 3 and 3^ represent another example when two branches are matched 
so that they form an inverted "seagull" . On the basis of these "elementary" solutions, various 
complex compound solutions can be constructed including periodic or chaotic stationary waves. 

The dashed line 1 in the figure corresponds to = (Q = 8). Another branch of the solution 
with the same value of £" = represents a solution of positive polarity which is unbounded both 
on ^ and y. For positive values of solutions of negative polarity become wider and of greater 
"amplitude" (see line 2). When E further increases and approaches Pmax, the solution becomes 
infinitely wide, but its minimum goes to —2. In the limiting case E = Pmax {Q = 1) two 
independent branches of the solution can be matched differently as shown by dashed-dotted 
lines 3 and 3^ in Fig. 9. The solution vanishes in this case when ^ = and goes to —2 when 
^ ioc; this situation is described by equation (13). 

For the negative E there are two families of solutions: negative one, corresponding to the 
left-hand side trajectories in Fig. 8, and positive one, corresponding to the right-hand side 
trajectories. When being negative, increases in absolute value {Q varies from — oc to 0_), 
solutions depart from the line 1 in Fig. 9 and gradually squeeze to the origin (see line 4 for 
instance). For the same values of negative positive solutions originated at infinity also 
gradually shrink and collapse in the origin (lines 6 and 5 demonstrate this tendency). 

Consider now the case iii) shown in Fig. 7b. The potential function in this case has only one 
well of a finite depth so that Pmin = Cf/{4:C2) at it; = — 2C2/C1, where C2 is negative now. There 
are no bounded solutions for negative w; they exist however for positive w and E varying in the 
range Pmin ^ E < 0. The finite value of the potential minimum corresponds to the equilibrium 
point of the centre type in the phase plane. There is also a family of closed trajectories for 
the above indicated range of E variation (see Fig. 10); these trajectories correspond to periodic 
solutions. 
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Figure 9. Various solutions described by equations (10), (12), (13) and (15). Compactons of negative 
polarity: line 1: Q = oo; lines 2 and 2' \ (5 = 2; lines 3 and 3^ (5 = 1: line 4: (5 = — 0.1. Compactons of 
positive polarity: line 5: (5 = —0.1; line 6: (5 = —0.25. 

■fw' 




Figure 10. Phase portrait of equations (3), (4) for the case 3c iii) (only those trajectories are shown 
which correspond to particle motion within the potential well in Fig. 7b). The dot at the center of closed 
lines indicates an equilibrium point corresponding to the potential minimum {E = —2.5 for the chosen 
set of parameters: Ci = 1, (^2 = —0.1); line 1: E = —2; line 2: E = —1.5; line 3: E = —1; line 4: 
E = -0.5. 

As usual, closed trajectories around the center (E > Pmm) correspond to quasi-sinusoidal 
solutions. Whereas other closed trajectories (E > Pmm) correspond to non-sinusoidal periodic 
waves with smooth crests and sharp narrow troughs. The larger is the value of E^ the longer is 
the wave period. The period tends to infinity when E ^ The analytical form of this family 
of solutions is described by the following equation: 



ay) = ±2V^ 



V4Q(Q - 1) - (y + 2Q)2 
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Figure 11. Various solutions described by equation (16). Line 1 (quasi-sinusoidal wave): Q = 1.01; 
line 2: Q = 1.5; line 3: Q = 2; line A: Q = 2.5. Dashed lines shows the equilibrium state y = —2. 



+ 2Q arctan 



y + 2Q 



V4Q(Q-l)-(y + 2Q)2 



(16) 



where ^ = xV^^{Ci / €2)^ , y = w{Ci/C2), Q = Cj/{4C2E). Solution (16) is shown in Fig. 11 
for different values of Q (note that the solution is negative in terms of y because C2 < 0). As 
follows from equation (16), y varies in the range: 



Q + VQiQ - 1) 



<y < -2 



Q-VQ{Q-^) 



whereas the dependence of wave period A on Q is: A{Q) = 87vQy/Q. The wave period varies 
from Stt to infinity when Q increases from unity to infinity. 

From the extreme values of y (see above indicated range of its variabihty) one can deduce 
the dependences of wave maximum and minimum on speed in the original variables. The corre- 
sponding formulae are: 



u 



max, mm 



iv)^v + 



Ci 



pV + p 



T2C21 



1 



2C2ipV + p) 



2C2{pV + p) 



+ 1 



(17) 



where the upper sign in front of the root corresponds to the wave maximum and lower sign - 
to the wave minimum. These dependences are plotted in Fig. 12 for V > —P/p in accordance 
with the chosen values of constants Ci — 1^ C2 — —0.1 and £" < 0. The asymptote V — —jS/p is 
shown in the figure by the vertical dashed line. As follows from equation (15), wave maximum 
cannot be less than the certain value, L^max, which occurs at some speed Vi shown in Fig. 12. 

For all possible values of wave maximum iXmax > ^^max, two values of wave speed are pos- 
sible, i.e. two periodic waves of the same maximum (but not minimum!) can propagate with 
different speeds. This is illustrated by the horizontal dashed line shown in Fig. 12 and drawn 
for iXmax = 2.5. In original variables quasi-sinusoidal waves exist when the speed is close to its 
limiting value l^max — ~p (^2 ~^ ^'^^^^ waves with greater speed. When V < l^max, 

the wave minimum and maximum are close to each other. Then, when the speed decreases, the 
gap between wave maximum and minimum gradually increases and goes to infinity when the 
speed approaches its minimum value Vmin = —f^/p- 
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Figure 12. Dependences of wave maximum (solid line) and minimum (dashed-dotted line) on speed in 
the original variables, equation (17), as follows from the solution (16). Dashed vertical line corresponds 
to V = —/3/p. The plot is generated for Ci = 1, C2 = —0.1 and p = /3 = 1. 



4 General case: p -\- q 

Consider now a more general case when the coefficients in equation (1) are such that p + g ^ 0. 
The basic equation (4) can be presented in the new variables 77 = (p + ^)x/6 and v = {p^q)w/6 
with the same constant of integration E = —{pV + /3)/2, but with new effective potential 
function 

P(v)-v-^-^. (18) 

The potential function is monotonic when Ci — C2 — 0, and there are no bounded solutions 
in this case. Bounded solutions may exist if at least one of these constants is nonzero. Below we 
present possible forms of the potential function and corresponding phase portraits of bounded 
solutions for various relationships between constants Ci and C2. Qualitatively all these cases 
are similar to those which have been described already in the previous section, therefore we 
omit the detailed analysis and do not present analytical solutions as they can be obtained 
straightforwardly and expressed in terms of elliptic functions. 

4a. If C2 — 0, the potential function represents a set of antisymmetric hyperbolas located 
either in the first and third quadrants when Ci = — 1, or in the second and fourth quadrants 
when Ci — 1; this is shown in Fig. 13. 

For the case of Ci = 1 only bounded solutions of a compact on type are possible for positive v. 
Such solutions correspond to the motion of particle c shown in the figure down to the potential 
well. This family of pulse- type solutions exist both for negative and positive E] all of them are 
bounded from the top with the maximum values depending on have zero minimum values 
and infinite derivatives when v — {). Corresponding phase plane is presented in Fig. 14a. 

For the case of Ci = — 1 there are two possibilities: i) there is a family of compacton-type 
solutions with < 0; they correspond to the motion of the particle h down to the potential well 
(particle motion to the left from the top of the "hill" corresponds to unbounded solutions). Pos- 
sible values of particle energy E vary for such motions from minus infinity to Pmax = — 2^/— Ci, 
where Pmax is the local maximum of the lower branch of the potential function (see Fig. 13). 
The phase portrait of such motions is shown in the left half of the phase plane in Fig. 14b. 

ii) Another possibility appears for the particle motion within the potential well shown in 
the first quadrant of Fig. 13 (see the particle a). Within this well all phase trajectories are 
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Figure 13. Potential function for the case p + C2 = and two values of Ci: Ci = 1 (solid line), 

and Ci = —1 (dashed line). Dots a, b and c illustrate possible motion of a point particle in the potential 
field. 




Figure 14. a) Phase plane corresponding to the potential function with Ci = 1 and various values of E. 
Line 1: E = —2; lines 2: E = —1; lines 3: E = 0; lines 4: E = 1; lines 5: E = 2. All trajectories 
in the left half-plane correspond to unbounded solutions, b) Phase plane corresponding to the potential 
function with Ci = —1 and various values of E. Line 1: E = —4; lines 2: E = —3; lines 3: E = —2; 
lines 4: E = -1; lines 5: E = 2.1; lines 6: E = 2.5; lines 7: E = 3. 

closed and corresponding solutions are bounded and periodical; they can be expressed in terms 
of elliptic functions. The phase portrait of such motions is shown in the right half of the phase 
plane in Fig. 14b. 

4b. If Ci = 0, but C2 ^ 0, the potential function also represents a set of antisymmetric 
hyperbolas located either in the third quadrant and right half-plane in Fig. 15 when C2 = 1, or 
in the first quadrant and left half-plane of that figure when C2 = — 1. 

For the case of C2 = 1 there are two possibilities: i) there is a family of compact on-type 
solutions with v < 0; they correspond to the motion of the particle b down to the potential well 
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Figure 15. Potential function for the case p + <7 7^ 0, Ci =0 and two values of C2: C2 = 1 (solid 
lines), and C2 = — 1 (dashed lines). Dots a, 6, c and d illustrate possible motion of a point particle in the 
potential field. 




a) 



b) 



Figure 16. Phase plane corresponding to the potential function (18) with Ci =0. a) C2 = 1 and various 
values of E. Lines I: E = -3; lines 2: E = -2; lines 3: £; = -1.89; lines A: E = -1; lines 5: £; = 0; 
lines 6: E = 1] lines 7: E = 2. b) C2 = —1 and various values of E. Line 1: E = —1; line 2: £^ = 0; 
line 3: E = 1] lines 4: E = 2; lines 5: E = 2.5] lines 6: ^ = 3. All trajectories in the left half-plane 
correspond to unbounded solutions. 



(particle motion to the left from the top of the "hill" corresponds to unbounded solutions). Possi- 
ble values of particle energy E vary for such motions from minus infinity to Pmax = 3(— (72/4)-'^/^, 
where Pmax is the local maximum of the left branch of the potential function (see Fig. 15). The 
phase portrait of such motions is shown in the left half-plane in Fig. 16a. 

ii) Another family of compacton-type solutions exist with v > 0; they correspond to the 
motion of the particle c down to the potential well. Possible values of particle energy E for such 
motions vary from minus to plus infinity. Corresponding phase plane is presented in the right 
half-plane in Fig. 16a. 
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a) 



b) 




Figure 17. Potential function for the case p -\- q 0. a) Supercritical case: Ci = —1 and two values 
of C2: C2 = 1 (solid lines), and C2 = — 1 (dashed lines); b) marginal case: Ci = —3 and the same two 
values of C2; c) subcritical case: Ci = —5 and the same two values of C2 (in the last case the horizontal 
and vertical scales are doubled). 



For the case of C2 = — 1 bounded solutions are smooth periodic waves which correspond to 
the particle oscillations in the potential well shown in the first quadrant in Fig. 15. Energy is 
positive for such motion and varies from Pmin — 3 (—(72/4)-'^/^, where Pmin is the local minimum 
of the right branch of the potential function (see Fig. 15) to infinity. Analytical solution for such 
waves can be also expressed in terms of cumbersome elliptic functions. Corresponding phase 
plane is presented in Fig. 16b. 

4c. Consider now the case when both Ci ^ and C2 ^ 0. The shape of the potential 
function is more complex in this case in general and depends on the relationship between the 
constants Ci and C2. The number and values of the potential extrema are determined by the 
number of real roots of the equation P^{v) = 0, where prime denotes the derivative on v. This 
condition yields (see equation (18)): 

+ Civ + 2C2 = 0. 

For real constants Ci and C2 this equation always has at least one real root. The real root 

2 /3 

is single when Ci > Cf^ = — ; its value is given by the expression 



V — 



(Ci/3)3 + C| 



C2 



1/3 



(Ci/3) 



(Ci/3)3 + C|-C2 



-1/3 



For the case Ci > Ci^, possible qualitative configurations of the potential function are shown 
in Fig. 17a for the particular choices of constants: Ci = — 1 and C2 — ±1. There is only one 
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Figure 18. Phase plane corresponding to the marginal case, Ci = Cf. a) Ci = —3, C2 = —1. Line 1: 
E = -3.05; line 2: E = -3; line 3: E = -2.9; line 4: E = -2; line 5: E = 3.8; line 6: E = 4; line 7: 
E = 4.5. All trajectories in the left half-plane correspond to unbounded solutions, b) Ci = —3, C2 = 1. 
Line 1: E = -5; lines 2: E = -4; lines 3: ^ = -3.75; lines 4: E = -3.5; line 5: E = 2.75; line 6: E = 3] 
line 7: E = 3.25; line 8: ^ = 3.5. 



local minimum at the right branch of the potential function for C2 = — 1 and a local maximum 
at the left branch of the potential function for C2 = 1. Almost the same configuration of the 
potential function occurs for the marginal case Ci = CJ^, as shown in Fig. 17b, however one 
more local extremum appears - on the left branch when C2 = — 1 and on the right branch when 
C2 = 1. In the case Ci < Cf^ the potential function is shown in Fig. 17c; there are three local 
extrema of the potential function for any value of C2 = ±1. 

The potential configuration in the supercritical case Ci > Cf^ qualitatively is similar to the 
case shown in Fig. 15, therefore the corresponding phase portraits are similar to those shown in 
Fig. 16. In the marginal case, Ci = Cf^, the potential configuration is also similar to those two 
cases mentioned above, however there are some peculiarities in the phase planes reflecting the 
appearance of embryos of new equilibrium points. Corresponding phase portraits are shown in 
Fig. 18. The embryos appear in the vicinity of = — 3 in Fig. 18a and in the vicinity of = 3 
in Fig. 18b. 

In the subcritical case Ci < C^^ the situation is different from the previous ones and should 
be considered separately. In the case of C2 = — 1, there are two potential wells, one of a finite 
depth on the left branch of function P(v) and another infinitely deep and wide well but bounded 
from the bottom on the right branch of function P{v) (see Fig. 17c). 

For the first potential well there is a family of closed trajectories in the phase plane corre- 
sponding to periodic solutions with the parameter E varying between the local minimum and 
maximum of the potential function; these solutions are described by elliptic functions. All closed 
trajectories are bounded by the loop of separatrix designated by symbol 3 in Fig. 19a. Trajecto- 
ries inside the separatrix loop next to center correspond to quasi-sinusoidal waves, and the loop 
of the separatrix corresponds to the solitary wave (soliton) which can be treated as the limiting 
case of periodic waves. The soliton shape is described by the following implicit formula: 



V 



±V2 



vi 



^/V2 -Vl + a/u2~ 



In 



(19) 
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a) 



b) 



Figure 19. Phase plane corresponding to the subcritical case Ci < . a) Ci = —5, C2 = —1. Line 1: 
E = -6; hnes 2: E = -5; hne 3: E = -4.25; hnes 4: E = -4; hne 5: E = 4.7; hne 6: E = 5; hne 7: 
E = 6; hne 8: E = 8. Ah trajectories in the left half-plane outside of the closed loop of separatrix 
correspond to unbounded solutions, b) Ci = — 5, C2 = 1. Line 1: E = —7; line 2: E = —5; lines 3: 
E = -4.657; lines 4: E = -4; lines 5: E = 4.3; line 6: E = 5] lines 7: E = 6.656; lines 8: E = 7. 





a) b) 
Figure 20. Soliton solutions on pedestals as described by equations (19). 



where vi^2 = -{CiT ^C'f - 3EC2) / E, {vi < V2) and E = PmUCi. where Pniax(Ci, C2) is 
the value of the potential local maximum shown in the left half-plane of Fig. 17c. Solution (19) 
is shown in Fig. 20a. 

In original variables function u describing soliton varies in the range 

V+^<u<V+ 



p + q p + q' 

thus, the soliton amplitude amounts 



^^^Q ^2-Vl ^ 12 ^Cl-3C2PmUCl.C2) ^ 
P + q p + q i^max(C'i,C2) 

The soliton velocity is 



(20) 



P 



(21) 
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Equations (20) and (21) allow one to obtain a direct relationship between the soliton's velocity 
and amplitude: 



For the second potential well located in the right half-plane of Fig. 17c, there is another 
family of closed trajectories in the phase plane corresponding to periodic solutions with the 
parameter E varying between the local minimum of the potential function and infinity; these 
trajectories are shown in the right half-plane of Fig. 19a. 

In the case of C2 = 1, there is a shallow potential well on the right branch of function P{v) and 
one infinitely deep well at the origin where the potential function is singular. For the shallow well 
there is a family of closed trajectories in the phase plane corresponding to periodic solutions with 
the parameter E varying between the local minimum and maximum of the potential function. 
All such trajectories are also bounded by the loop of separatrix designated by symbol 7 in 
Fig. 19b. The loop of separatrix corresponds to the solitary wave whose shape is described by the 
same implicit formula (19), but with different values of constants C2, E and Pmax(C'i, C2), 
where Pmax(C'i, C2) is the value of the potential local maximum shown in the right half-plane of 
Fig. 17c. This solution is shown in Fig. 20b. All above relationships between soliton amplitude 
and velocity, as well as between soliton amplitude or velocity and constants Ci and C2 remain 
the same as above. 

For the infinitely deep well at the origin there are two families of compactons with nonpositive 
and nonnegative values; the phase plane for them is similar to that shown in Fig. 8 and solutions 
are similar to those shown in Fig. 9. The entire phase portrait of the system in the case of C2 = 1 
is shown in Fig. 19b. Phase trajectories corresponding to positive compactons are not shown in 
detail in that figure because they are too close to each other and are in the narrow gap between 
the axis and external two unclosed branches of the separatix 7 (only one such trajectory, 
line 5, is shown in Fig. 19b; all other trajectories are similar). 

5 Conclusion 

As was shown in the paper, the extended reduced Ostrovsky equation (1) possesses periodic 
and solitary type solutions in general. There is a variety of solitary-wave solutions including 
compactons with infinite derivatives at the edges, smooth solitons, and periodic waves. All 
compactons, however, actually are of the compound- type solutions, i.e., they consist of two 
or more non-smooth branches. Among periodic waves depending on the equation parameters, 
there are also both smooth solutions and compound-type solutions which consist of periodic 
sequences of non-smooth branches (see, e.g.. Fig. 5b). Moreover, using compacton solutions as 
the elementary blocks, one can construct very complex compound solutions including stochastic 
stationary waves. 

The approach used in this paper and based on the qualitative theory of dynamical systems 
is free from the limitations of paper [1] and allows us to present a complete classification of all 
possible solutions of stationary exROE. In particular, solutions were obtained and analyzed in 
details for the case p -\- q = that was out of consideration in paper [1]. Another "prohibi- 
ted" combination of parameters, qV — 7^ 0, that was also out of consideration in paper [1], 
does not even appear in our study. The approach exploited in the present paper is based on 
a vivid mechanical analogy between a particle moving in a special potential field and considered 
stationary exROE. This approach allows one to observe qualitatively an entire family of all 
possible solutions even without construction of particular exact solution. A similar approach has 
been exploited recently in application to the reduced Ostrovsky equation [4, 3] and exROE [3], 
although in the last case, the complete solution classification was not considered. 
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